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ABSTRACT 

We demonstrate how the Fundamental Manifold (FM) can be used to cross-calibrate distance es- 
timators even when those "standard candles" are not found in the same galaxy. Such an approach 
greatly increases the number of distance measurements that can be utilized to check for systematic 
distance errors and the types of estimators that can be compared. Here we compare distances ob- 
tained using SN la, Cepheids, surface brightness fluctuations, the luminosity of the tip of the red 
giant branch, circumnuclear masers, eclipsing binaries, RR Lyrae stars, and the planetary nebulae 
luminosity functions. We find no significant discrepancies (differences are < 2a) between distance 
methods, although differences at the ~ 10% level cannot yet be ruled out. The potential exists for 
significant refinement because the data used here are heterogeneous B-band magnitudes that will 
soon be supplanted by homogeneous, near-IR magnitudes. We illustrate the use of FM distances to 
1) revisit the question of the metallicity sensitivity of various estimators, confirming the dependence 
of SN la distances on host galaxy metallicity, and 2) provide an alternative calibration of Ho that 
replaces the classical ladder approach in the use of extragalactic distance estimators with one that 
utilizes data over a wide range of distances simultaneously. 



1. INTRODUCTION 

Distance determinations to galaxies are key to many 
areas of study in extragalactic astronomy, but remain 
problematic. The potential for systematic errors lingers 
above much of the work, particularly at the accuracy 
level of a percent or less cu rrently desired (see review by 
iFreedman fc Madord l2010h . These concerns can be ad- 
dressed through cross calibration among techniques that 
have independent challenges. For example, one estima- 
tor might have a hidden dependence on metallicity, an- 
other might be more susceptible to internal extinction. 
One state-of-the-art recent study compares distances ob- 
tained using SN la and Cephe id variable stars in the 
same galaxies (jRiess et al.ll201ll ) to claim a distance cal- 
ibration precision of a few percent. 

Such cases, where one has distance measurements from 
a variety of methods for a specific galaxy, are rare for 
two reasons. First, because distance estimators are of- 
ten related to distinct stellar populations, they are not 
often found within the same stellar system. The classic 
example is the dichotomy between t he Galactic e nviron- 
ments of RR Lyrae and Cepheids flBaadd [1944D . Sec- 
ond, and currently most relevant, is that several factors 
conspire to limit estimators to distinct distance regimes. 
Cepheid observations, like those of RR Lyrae, eclipsing 
binaries and the tip of the red giant branch (TRGB), re- 
quire observations that resolve the stellar populations of 
the system and are therefore limited to relatively nearby 
galaxies, even when using the high angular resolution 
provided by the Hubble Space Telescope. On the other 
hand, SN la, due to their rarity, tend to be found in more 
distant galaxies simply because of the gain in search vol- 
ume with increasing distance. These effects make it quite 
difficult to compare the r esults among certa in classes of 
distance estimators. The lRiess et al.l ([20111 ) study cited 
above, which is among the most extensive and complete, 
is based on only eight galaxies with both Cepheid and 



SN la measurements. 

One solution to this limitation relies on identify- 
ing a distance estimator that can be applied to all 
galaxies, near and far, early or late type, giant or 
dwarf. Even if such an estimator has a large intrin- 
sic distance uncertainty for any individual galaxy, and 
so is not ultimately superior to the other estimators, 
the gains reaped by being able to compare the rela- 
tive precision and accuracy of the preferred estimators 
over many tens, if not hundreds of galaxies, is criti- 
cal. Distance det erminations from scal ing relations, such 
as Tull y-Fisher (iTulry fc Fisher] 119771) or Fundamenta l 
Plane ([Dressier et al.l Il987t iDiorgovski fc Davisl 119871) . 
span the range of distances probed by many methods 
and so provide a fiducial against which different distance 
estimators can be compared. Their shortcoming is that 
they are each applicable to different, narrow classes of 
galaxies, limited in type and luminosity, and therefore 
do not resolve the entire problem. 

In a series of papers, we have presented a new scal- 
ing relation that is applicable to all ga laxies, regard- 
less of lu minos i ty or morphological class (Zaritsk y et al.1 
I2006albl 12001 |20T1 . This relation, termed the Fun- 
damental Manifold (FM) in reference to its most direct 
antecedent, the Fundamental Plane, presents an oppor- 
tunity to uncover hidden systematic differences among 
the current set of distance estimators. Here we use exist- 
ing data, distances from the NED 1-D database and pho- 
tometry a nd kinematics comp iled within the Hyperleda 
database (|Paturel et al.l [20031) . to test this basic idea. 
Neither the data, due to their inhomogeneity, nor the 
sample, due to the intractable selection criteria, are op- 
timal for precise cross-calibration of distance estimators 
and eventual application to related problems, such as 
the measurement of Ho . Therefore, the results presented 
here should be viewed as only illustrative of what an op- 
timized treatment using the FM could yield. In §2 we 
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Fig. 1. — Distribution of distances as a function of distance- 
measuring technique among NED-ID sample (from top, super- 
nova la peak luminosity corrected for decay rate, planetary neb- 
ular luminosity function, surface brightness fluctuations, Cepheid 
period-luminosity relation, the luminosity of the tip of the red gi- 
ant branch, circumnuclear masers, eclipsing binaries, and the RR 
Lyrae period-luminosity relation). Y-axis arbitrarily normalized 
to match scales among subsamples. Although some distance es- 
timators span the same range of distances, the differences among 
distance ranges can be as large as two orders of magnitude. 

discuss our selection of the sample, in §3 we present the 
results of this comparison among distance estimators, il- 
lustrate how the FM can be used to trace potential inter- 
nal systematic dependencies, and reverse the calibration 
and use the FM to obtain distances and estimate Hq. In 
§4 we summarize our results. 

2. THE DATA 

2.1. Distances 

We utilize the NED 1-D database, which is a compi- 
lation of published distances from a wide ranging set of 
methods. Without prejudice for or against any of the 
entries in the database, we calculate the weighted aver- 
age (weighting inversely by the quoted uncertainty and 
neglecting the measurement if no uncertainty is given) 
for each of the distance estimators, for each galaxy. We 
require at least two measurements using a particular 
method to calculate a mean distance measurement for 
a given galaxy. There are 1153 galaxies for which we 
retrieve a distance using at least one of the methods we 
are comparing. In Figure [T] we show how distance estima- 
tors are often effectively constrained to distinct distance 
ranges. More critically, the distribution of objects with 
SN la distance measurements overlaps very little with the 
other methods, limiting the number of systems that can 
be used to directly check for potential systematic errors 
in the cosmologically-critical la distances. 

2.2. Photometry and Kinematics 

To obtain the structural parameters necessary to 
place a g alaxy on the FM, we search the Hyperleda 
database (|Paturel et al.ll2003l ). which provides morpholo- 
gies, kinematic measurements, magnitudes and struc- 
tural parameters, the latter two based on the work 



Fig. 2. — Distribution of morphological types as a function of dis- 
tance measuring technique among NED-ID sample. Y-axis arbi- 
trarily normalized to match scales among subsamples. Differences 
among the types of galaxies accessible to various methods are clear. 

of iPrugniel &: Heraudeaul jT998), for those galaxies for 
which we obtained distances from NED 1-D. For the 
H I rotation measurements, we adopt the homogenized 
values provided in Hyperleda, which include an incli- 
nation correction. We calculate weighted averages of 
both the stellar rotation and velocity dispersion mea- 
surements where available, and only accept values only 
if there are at least two measurements of either quantity, 
and if either the rotation or dispersion is larger than 10 
km sec -1 , to avoid the systems with the largest frac- 
tional uncertainties. For the morphologies, we average 
T- Types using uniform weights and require a minimum 
of only one measurement. The photometric parameters 
are taken from the archival B-band data, which is the 
most prevalent. There are various reasons (star forma- 
tion, extinction) to use redder-band photometry, and we 
will discuss results using the available I -band data. In the 
future, we will use the S 4 G database (jSheth et al.ll2010l) 
to obtain photometric parameters at 3.6/xm, which ap- 
pear to significantly reduce t he scatter in scaling rel ations 
(jFreedman fc Madordl2mTl : iFreedman. et alJ[20ll . The 
distribution of morphologies, for the subsample of galax- 
ies for which morphologies are catalogued, is shown in 
Figure [3] Again, we find that certain distance estima- 
tors are limited in their coverage, making comparison on 
an object-by-object basis difficult. There are 343 galax- 
ies with both distance measurements and the necessary 
data to place them on the FM. 

For a subset of the data, we also o btain an optical color, 
(B — V)o, from the RC3 catalog (|de Vaucouleurs et al.l 
|1991|) to investigate the role of color, and, by inference, of 
differences in the stellar mass-to- light ratio, on a galaxy's 
deviation from the FM. We use observed magnitudes 
rather than k-corrected rest-frame magnitudes, because, 
despite the large range in distance of the NED 1-D sam- 
ple, those galaxies with all of the necessary data lie at 
cz < 6000 km sec -1 . One evident avenue for progress is 
obtaining these data for a larger fraction of the galax- 
ies with SN la distances, and thereby also extending the 
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sample to greater distances. 

Our compilation of the available data is presented in 
Table [T] (complete Table available electronically). 

3. THE FUNDAMENTAL MANIFOLD AND TESTS OF 
DISTANCE ESTIMATORS 

3.1. Overview of the Fundamental Manifold 

The FM has been de scribed in detail elsewhere 
(jZaritskv et alj 120081 1201 ID . but, in summary, it is de- 
rived from the virial theorem with the additional infor- 
mation that the mass-to-light ratio within r e , T e , can 
be either 1) measured independently (e.g., via strong 
gravitational lensing) or 2) estimated using an empiri- 
cal function of the internal kinematics, V, and the sur- 
face brightness within r e , I e . In the latter case, we have 
found that the structural differences (spatial and kine- 
matic) among galaxies, other than those captured in the 
variations in T e , are insignificant given the observational 
errors (~ 0.1 dex). The FM is a conceptually simpler 
scaling relation than either Tully-Fisher (TF) or Funda- 
mental Plane (FP) in that it can be expressed indepen- 
dent of wavelength and has only a zero point term to 
calibrate (rather than slope in TF or "tilt" in FP): 

log r e = 2 log V - log I e - log T e - C, (1) 

where V is expressed as the combination of pressure 
(a) and rotational support (v r ), V 2 = a 2 + v 2 /a {a 
is expected to lie be tween 2 and 3; We iner et al.1 12005 
iZaritskvet aT1l20rl l and C is a constant that is ef- 
fectively the zero point of the FM distance estima- 
tor. With measurements (or estimates) of all terms 
on the right hand side of the equation, one solves for 
r ei which, in concert with the angular measurement 
of r e , provides the distance. Because we generally 
lack independent estimates of T e , we use the empiri- 
cal relationship for T e (V,I e ) derived from existing data 
(jZaritskv et al.l 120111) . Use of the empirical relation- 
ship reintroduces wavelength dependences and additional 
complexity (analogous to FP "tilt" ) , but the FM still re- 
tains the advantage over previous scaling relations in that 
it is applicable to all galaxies. 

As with any empirical relationship, the relationship 
changes slightly depending on the data used to define 
it, and uncertainties in the fitting function ultimately 
translate to systematic uncertainties in the distances es- 
timated using the FM. Consequently, comparing FM dis- 
tances to those obtained using other distance estimators 
provides a test of both the FM and the independent dis- 
tance estimators. Comparisons of the FM and multiple 
independent methods will allow us to distinguish between 
systematic problems in the FM versus problems with 
the inditvidual independent estimators. Here we choose 
to define the fitting functio n using only the sph eroidal 
galaxies in the original data (Zaritskv et al. 2008?), which 
are expected to have minimal stellar mass-to-light ra- 
tio variations. We exclude globular clusters and ultra- 
compact dwarfs, wh ose nature is still controversial (see 
iZaritsky et al.ll2011l ) and that extend the FM beyond the 
r e range needed for this work. Although changes in the 
fitting function for T e will affect the results, the fitting 
function is most constrained in exactly the parameter 
range most relevant here — that of disk and spheroid 
galaxies with luminosities around L». The functional 
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Fig. 3. — The relationship between color and v r /a for galaxies 
in our sample that have published measurements for stellar ve- 
locity dispersion, <r, and rotation, v r . We select the subsample 
with B — V > 0.85 and v r /cr > 0.2, separated with dashed lines, 
for subsequent analysis because it is large and has no measurable 
relationship between color and v r /a (see text for an expanded dis- 
cussion of this selection). 

form we use is 

log T e = 2.12 - 0.28 log V - 0.82 log J e + 0.15 log 2 J^l) 

+0.26 log 2 !/ - 0.09 log V\og I e . (2) 

This equation provides an estimate of T e in the I-band. 
There are two related, but distinct, questions that will 
be addressed in the following section: how much does 
the use of B-band magnitudes increase the scatter (there 
would be no increase if all galaxies had the same B — I 
color) and how much scatter is introduced in the /-band 
T e by differences in stellar populations? 

To test each distance method, we use the distances 
provided by NED 1-D to calculate the physical value of 
r e and then place the galaxy on the FM. In our original 
work, from which we take our fitting function for T e , we 
used simple Hubble flow distances with an adopted Hub- 
ble parameter. Any difference between the true Hubble 
parameter and that adopted previously will affect the cal- 
ibration of T e and C, manifesting itself as a zero point 
shift (i.e., as a change in C). We therefore have freedom 
to adjust the zero point of the FM as needed (see below 
for our treatment of this) and so our discussion here fo- 
cuses on relative differences between distance estimators 
rather than on an absolute calibration. 

3.2. Morphology Related Differences 

One consideration before proceeding with our compar- 
ison of distance estimators and the FM is that of the 
appropriate choice of a in the definition of V . In previ- 
ous work we adopted a = 2, as also don e by others (e.g., 
IBurstein et al.l [19971 : IKassin et al.ll2007D . However, rea- 
sonable model s suggest that value s between 2 and 3 are 
plausible (see iWeiner et al.1 12006D . Because v r /a maps 
onto morphological type, and morphological type maps 
onto color, which is presumably related to the stellar 
mass-to-light ratio, there is the potential to misidentify 
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TABLE 1 

The Sample 



Name 




(B-V)o 


He 


2r e 
["] 


[km/s] 


Vr,* 

[km /s] 


a 

[km / s] 


Dc 
[Mpc] 


D T 
[Mpc] 


D Ja 

[Mpc] 


DsBF 

[Mpc] 


DpjV e 

[Mpc] 


Djjj? 
[Mpc] 


Dm 
[Mpc] 


D S4 
[Mpc] 


cz 
[km/s] 


NGC7814 


11.38 


0.83 


22.23 


118.1 


230.7 


0.0 


172.9 








12.70 










1047 


NGC0055 


8.68 


0.54 


22.18 


400.1 


58.9 


0.0 


0.0 


1.92 


2.10 














142 


NGC0147 


10.65 


0.78 


24.02 


376.8 


0.0 


0.0 


26.4 




0.64 








0.63 






-176 


NGC0150 


11.90 


0.50 


21.40 


63.4 


164.9 


0.0 


0.0 


















1566 


NGC0185 


10.38 


0.73 


22.11 


177.1 


0.0 


0.0 


23.3 




0.61 














-201 


NGC0205 


9.00 


0.82 


21.82 


292.5 


0.0 


22.0 


33.5 




0.80 














-258 


NGC0221 


9.19 


0.88 


18.75 


65.2 


0.0 


37.8 


80.9 




0.69 




0.77 




0.82 






-203 


NGC0224 


4.88 


0.68 


19.89 


802.0 


0.0 


75.1 


183.2 


0.80 


0.79 




0.75 


0.72 


0.77 




0.75 


-294 


NGC0247 


9.42 


0.54 


23.37 


492.2 


249.8 


0.0 


0.0 


3.76 


3.60 














142 


NGC0300 


8.51 


0.58 


22.87 


594.5 


85.1 


0.0 


0.0 


1.99 


1.96 
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a relationship between the FM residual (the deviation 
of a galaxy from the FM) and a, with one that is truly 
due to color (or vice- versa). Without the proper choice 
of a and some means of avoiding differences in stellar 
mass-to-light ratios beyond those included in the empiri- 
cal calibration of T e , there will be systematic differences 
when comparing distances to early and late type galax- 
ies. We will mitigate this difficulty by 1) proceeding to 
optimize our choice of a, 2) removing any residual cor- 
relation between FM residuals and color, and 3) provid- 
ing comparisons among galaxies of similar morphological 
type. 

There is indeed a strong correlation between (B — V)q 
and v r ja (a Spearman rank correlation test for our 
sample puts the probability of it arising at random at 
4 x 10~ 5 ). As such, for the purpose of constraining a we 
work with a subsample for which v r /a is independent of 
color (Figure H; {B - V) Q > 0.85 and v r /a > 0.2 (the 
probability of the correlation arising randomly for this 
sample is 0.79). This particular part of the analysis, the 
derivation of a, is restricted to early type galaxies be- 
cause a is generally not available for disk galaxies. As 
a result, the galaxies in Figure |3] are naturally restricted 
to fairly red colors. If a hidden correlation remains, we 
expect some of that to influence the best fit value of a, 
again because v/a is expected to correlate with color, 
and eventually to lead potentially to differences between 
FM residuals and morphological type. 

For this subsample, we calculate the root mean square 
(rms) residual about the FM (allowing the zero point to 
float) and identify the value of a that minimizes the rms 
(a = 2.68; see Figure EJ. Reassuringly, this value lies 
within the theoretically plausible range and we adopt 
this value of a for the remainder of this study. How- 
ever, this calculation is unfortunately quite sensitive to 
the selected subsample. For example, if we select galax- 
ies that have (B — V)q > 0.8 (instead of 0.85), then the 
minimum rms residual occurs at a — 1.81 (although a 
correlation between color and v r /a still exists at > la 
confidence for this sample and so this sample is not pre- 
ferred). The robustness of the a determination needs to 
be improved both by obtaining a larger sample of galax- 
ies and by performing this analysis with photometry in 
a filter band that is less sensitive to morphological type 
(and, by relation, to color). The effect of varying a is 
principally to alter the relative offset between early and 
late type galaxies relative to the FM (see §4), but does 




Fig. 4. — Residual about the FM as a function of adopted a. The 
minimum value at 2.68 is chosen for the remainder of the study. 



not affect the conclusions presented in §5. 

Any remaining systematic behavior of the FM residu- 
als will point to new phenomenon, including systematic 
errors in distance estimators. We begin by examining the 
dependence of FM residual to color (Figure [5]). The data 
show a trend of residual with color. It is not surpris- 
ing that the early-types have near zero mean residual, 
because the FM relationship we use was defined using 
spheroidal galaxies (f rom a different sample based on a 
range of sources, see iZaritskv et al.1 (|2006al) . but early- 
types nonetheless). Neither the sense nor the magnitude 
of the trend is necessarily straightforward to explain. 
There is an expectation of a trend with color because the 
stellar mass-to-light ratio will change with color. How- 
ever, other factors change as well. For example, the gas 
mass fraction of a galaxy correlates with color and that 
mass is unaccounted for in the empirical derivation of T e 
because we used only spheroidal systems. 

We use the mean relationship, shown as the line, to 
correct all of our data (and so subsequently restrict our- 
selves to systems with available color measurements). If 
a particular estimator (including the FM) has a system- 
atic error that depends on galaxy type (color), it will 
contaminate this plot. As with many of the issues raised, 
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we expect this problem to be mitigated using redder pho- 
tometry, but proceed nonetheless with the data currently 
available. Recall that our expression for T e provides an 
estimate of the total /—band mass-to-light ratio within 
r e , even though we are using B— band magnitudes. This 
means that there are two potential sources of scatter in- 
troduced by stellar population differences. First, there 
is the possibility of large variations in B — I within the 
sample that then introduce scatter because we implicitly 
assume a constant B — I color (solar) when converting 
B-band luminosity into the corresponding mass in solar 
masses. Second, scatter can arise from stellar population 
differences and the associated variations in the stellar 
mass-to-light ratio. 

We show, in Figure [6l the distribution of B — I colors 
for our galaxy sample, and note that a large fraction of 
the galaxies are fairly well constrained to colors that vary 
by about 0.3 mag (a constant difference from solar color 
is not a cause for concern because we eventually renor- 
malize the entire relation using the SBF distances) . More 
importantly, this range of colors produces small scatter 
(0.046 dex) in T e between the values calculated using the 
/3-band magnitudes (and the color correction described 
above) and the /-band magnitudes (with a color correc- 
tion derived for the /-band). We conclude that varia- 
tions in B — I colors, beyond those that we remove using 
a mean trend of FM residual with (B — V)o, are a minor 
contributor to the scatter we observe in the FM (> 0.1 
dex for these data). 

Our (B — V)o color correction is intended to empir- 
ically remove such a correlation, although as hinted at 
previously, this correction may be removing other fac- 
tors as well that work in parallel to stellar population 
differences. In Figure [5] we plot the expe cted variation in 
stellar mass-to-light ratio vs. B — V from lBell fc de Jong! 
(120011 ). selecting one particular model (closed box, scaled 
Salpeter IMF), demonstrating that the effect we are re- 
moving is indeed of the same scale as that expected from 
stellar population variations. 

Finally, we estimate the degree to which going to even 
redder passbands might reduce the scatter, when one 
does not have a color measurement. For the same stellar 
population model as in Figure[SJ the slope in the logarith- 
mic relation between B — V and the stellar mass-to-light 
ratio at K is about half that in /. Therefore, if the scatter 
in the FM is dominated by stellar population differences, 
and we do not have colors, we still expect to reduce the 
scatter by two simply by observing at K. Further gains 
are p ossible if extinction is important (jFreedman. et ail 
1201 If ). If other causes, not related to the photometric 
passband, contribute comparably to the scatter, or if we 
can partially account for stellar population differences, 
then we expect less than a factor of two improvement 
by going into the infrared. For example, the scatter in 
our /3-band FM is the same as in our /-band FM (~ 0.1 
dex), because we have accounted for the stellar popula- 
tion differences at B. 

The key to interpreting our results is to determine 
whether multiple estimators show the same systematic 
behavior in FM residuals, in which case the problem is 
most likely to lie with the FM or the various corrections 
we have implemented, or if the FM residuals are pecu- 
liar for specific estimators, in which case the fault lies 
with that particular estimator. On the other hand, if 
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Fig. 5. — FM residual for individual galaxies as a function of 
color. Pressure supported galaxies (those with measured cr) are 
plotted as dots and represent the results using the distances from 
surface brightness fluctuations. For the rotationally supported 
galaxies, plotted as filled circles, we show the results using distances 
derived from Cepheid measurements, tip of the red giant branch, 
SNe la, planetary nebulae, and surface brightness fluctuations. The 
solid line represents the best fit and is used for the color-dependent 
correction. The dashed line represents the expected color depen- 
dence, normalized to cross our fit at (B — V)o = 0.7, derived from 
a synthetic stellar population (see text for details). 
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Fig. 6. — The effect of B — I color variations on the scatter in 
the FM. Left panel shows the distribution of B — I for the galaxies 
in our sample. The right panel shows the resulting differences in 
estimates of T e using the B and / band magnitudes, after removing 
for B — V color dependence. This comparison, which is unphysical, 
highlights the scatter introduces solely by assuming a uniform B — I 
color for all galaxies, but not assuming that the stellar populations 
are all the same. We conclude that the use of i?-band magnitudes 
is not a major source of scatter in the FM, if one does correct for 
differences in stellar mass-to-light ratios using colors. 
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Fig. 7. — FM using measured distances. We compare the FM's using different distance estimators and divide the sampie into pressure 
supported (upper panels) and rotationally supported (lower panels). Furthermore, we color and shape code according to galaxy morphology 
(blue squares for late type, red circles for early type). The x axis is log r e in kpc, and the y-axis is 21ogV — log/ e — logT e — C, where V 
is in km sec -1 , I e is in Lir,/sq. pc, and T e is in solar units. C is obtained by calibrating to the sample of surface brightness fluctuation 
distances (SBF) for the pressure supported galaxies. The line is the 1:1 expectation, not a fit to the data. Individual uncertainties are not 
plotted, but the scatter is used to estimate the uncertainty in the mean offsets from the 1:1 line. 



none of the estimators exhibit systematic problems rela- 
tive to the FM, then we will conclude that the estimators 
and the FM are accurate distance estimators to at least 
within the quoted uncertainties. 

3.3. Comparing Distance Estimators and the FM 

Our aim is to test a variety of popular distance esti- 
mators to determine which, if any, may lead to distance 
estimates with systematic errors. Using the FM as a 
reference, we now compare estimators that we could not 
previously test against each other. To do this, we plot the 
galaxies on the FM using the mean distances drawn from 
the NED 1-D data (Figure [7]). Normalization or slope er- 
rors for any particular estimator will suggest that there 
is a problem with that estimator. On the other hand, if 
the normalization or slope errors appear endemic, then 



the problem will lie with the FM. 

For each distance estimator with more than 10 galaxies 
with all the necessary data, we present the results, show- 
ing the pressure supported systems, which have velocity 
dispersions and sometimes also have stellar rotation val- 
ues, v r (upper panels in the Figure), and the purely rota- 
tionally supported galaxies, which have no quoted stellar 
velocity dispersion and for which we adopt the H I width 
as a measure of v r (lower panels). In general, if there is a 
stellar velocity dispersion measurement, the galaxy is a 
pressure supported systems (v r ja < 1). Within each 
panel we also color-code based on morphology, divid- 
ing the early and late type galaxies at a T-Type of 1. 
We have set the zero point of the FM, C, using the re- 
sults from the pressure-supported, SBF sample, which is 
our largest subsample of galaxies and also predominantly 
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Fig. 8. — FM using measured distances. Same as FigureQ except 
we include other distance estimators (maser (blue square), eclipsing 
binaries (red triangle, plotted large to enhance visibility), and RR 
Lyrae (blue circles)). The value of C is as in Figure [71 

consists of early-type galaxies, which will be less suscep- 
tible to extinction and stellar population variations. 

A cursory examination of the panels shows that there 
are no disastrous problems with any of the distance es- 
timators as applied to any of the galaxy subsamples, 
even though some perform better than others (either in 
terms of zero point or scatter). There are a few individ- 
ual galaxy outliers, although it is often the same galaxy 
that is an outlier in multiple panels because distances 
are available from multiple estimators, suggesting that 
fault lies not with the distance estimate but rather with 
one of the other parameters that enters the FM. NGC 
4704, which has an unusually large residual in the SN 
la panel, is removed from subsequent analysis of the SN 
la distance estimator. In Figure [8] we show the results 
for galaxies that have distances measured using eclipsing 
binaries, RR Lyrae, and masers. 

In Figure IH1 we compare the mean residual and the un- 
certainty in the mean between the various distance meth- 
ods, for pressure and rotationally supported galaxies (the 
values are given in Table [2]) . The agreement is strikingly 
good, with the possible exception of results obtained us- 
ing the "other" methods, although there the uncertainty 
is large because of the small number of galaxies in that 
category. None of the differences are statistically signifi- 
cant, but the magnitude of differences is in the range that 
is potentially important for distance estimators (~ few to 
tens of %). For example, the difference in zero points ob- 
tained using the TRGB and SN la methods and pressure 
supported galaxies is 0.062, which corresponds to an in- 
ferred distance difference of 15%. Unfortunately, for the 
purpose of identifying if there are any real discrepancies 
at the target level of a few percent, the error bars are still 
somewhat too large due to the large amount of scatter 
in the FM (arising from the heterogeneous B-band data) 
and small sample size (as many galaxies lack the full set 
of measurements needed to calculate the FM). 

The uncertainties in our comparison of distance esti- 
mators are currently dominated by the number of galax- 



FlG. 9. — Mean offsets from FM using distances derived from dif- 
ferent estimators. Red circles denote the results for pressure sup- 
ported galaxies, while blue squares those for rotation supported 
galaxies. "Other" includes results from eclipsing binaries, RR 
Lyrae and masers for both rotation and pressure supported galax- 
ies. Error bars represent the dispersion in the mean value. The 
distances are referenced by definition to the surface brightness fluc- 
tuation (SBF) calibration for pressure supported galaxies. 



TABLE 2 

Mean Offsets from FM 



Distance Estimator Pressure Supported Rotationally Supported 



Cepheids -0.002 ± 0.134 -0.025 ± 0.023 

SBF 0.000 ± 0.009 -0.004 ± 0.039 

SN la -0.062 ±0.027 -0.020 ± 0.027 

TRGB 0.043 ± 0.068 0.003 ± 0.048 

PNe 0.014 ± 0.034 -0.049 ± 0.085 



ies in each class and the intrinsic scatter in the FM. 
The small sample sizes can be addressed relatively easily. 
Many galaxies with extant distance measurements, par- 
ticularly those with SN la distances, lack the additional 
measurements (typically V) necessary to place them on 
the FM. These can be obtained with ground-based ob- 
serving. The FM scatter can be improved by going to 
redder pass-bands. If all of the data had uncertainties 
similar to the pressure-supported, SBF measurements, 
we would easily confirm or refute the level of discrep- 
ancy among distance estimators found here. 

3.4. Testing the Metallicity Dependence of Distance 
Estimators 

The FM is not only of use in uncovering systematic 
differences among distance estimators, but is also po- 
tentially useful in uncovering systematic errors in any 
particular distance estimators within a set of galaxies. 
One long-running concern is the role chemical abun- 
dance might play in affecting distance estimat ors (see 
iFreedman fc Madord[201ll iGallagher et al.ll2008l for dis- 
cussions of the effects on Cepheids and SNe la, respec- 
tively). Given the correlation bet ween host galaxy gas - 
phase metallicity and luminosity (Zarit skv et al.|[l994D . 



Zaritsky et al. 



Cepheids 



1 
log r. [kpc] 



Fig. 10. — Residuals from the FM as a function of galaxy absolute 
magnitude for different distance estimators. As in other Figures, 
late type galaxies plotted as blue squares and early type galaxies 
as red circles. Systematic deviations, if present, suggest a possible 
abundance dependence in the respective distance estimator due 
to the correlation between galaxy luminosity and mean chemical 
abundance. 



Fig. 11. — The Fundamental Manifold relation for galaxies with 
distances determined with at least two independent methods. We 
adopt the average of the distance determinations. Pressure sup- 
ported systems are plotted as red circles and rotationally supported 
systems as blue squares. When we have both measurements of the 
HI circular velocity and stellar kinematics, we place the galaxy 
on the FM using each of these and connect the points. These 
connections illustrate the contribution of the uncertainties in the 
kinematics to the scatter. 



TABLE 3 

Magnitude Dependent Offsets from FM 



Distance Estimator 



All Data 



Faint Tail Trimmed 



Cepheids 
SBF 
SN la 
TRGB 
PNe 



-0.010 ±0.011 
0.033 ± 0.008 
0.013 ± 0.032 
0.007 ±0.015 

-0.002 ± 0.029 



0.023 ± 0.022 
0.037 ± 0.008 
0.093 ± 0.045 
-0.043 ± 0.042 
0.053 ± 0.035 



we use our sample to examine whether there are potential 
dependencies on metallicity by measuring correlations 
between FM residuals and luminosity (a correlation could 
also implicate another parameter that correlates with lu- 
minosity, e.g., mean stellar age). In Figure [TU] we show 
the FM residuals as a function of galaxy absolute magni- 
tude for our five primary estimators. Visually, two of the 
panels appear to show a trend of increasing FM residual 
with decreasing luminosity (SBF and SN la). However, 
this similarity among panels overstates the concordance 
because some of the galaxies, particularly at the bright 
end, appear in multiple panels. In Table [3l we present 
measurements of the relationship between the FM residu- 
als and absolute magnitude for each distance estimator. 
We list both the fit obtained using all of the galaxies 
and after trimming the poorly populated faint end of 
the galaxy distribution (either at Mb = —17 or —18). 
The only two significant correlations are provided by the 
SBF and SN la data (their rank correlation probabilities 
of arising at random are both < 0.3%, if we confine the 
SNe sample to galaxies with Mb < —18). 

If we attribute the slopes to a metallicity dependence 
in the estimator, we can quantify the effect and compare 
to previous determinations. Our best fit for the SN la 
in galaxies with Mb < —18 implies a distance offset of 



- 41 ± 0.11 dex/ [ Q/Hl (for the 4.5 mag/[0/H] relation 
of iZaritsky et al.l (119941) V wh ich agrees with the se nse 
found bv lHowell et all ((20091) and IKellv et all pOlOh in 
that one would underestimate the distance to SNe in 
more metal rich galaxies, because SNe are intrinsically 
more luminous after light curve correction. It is diffi- 
cult to quantitat ively compare our results to those of 
IKellv et all (f2QToh for' various reasons, including their use 
of stellar masses vs. our use of luminosities and our cal- 
ibration to gas-phase abundances. However, they find 
distance modulus variations that correspond to about 
0.15 over the 1 order of magnitude in mass. If we as- 
sume that the order of magnitude in mass correspond to 
the same magnitude range in luminosity, then we would 
predict a distance effect of 0.23 dex over the correspond- 
ing range, which in turn results in a change of 1.15 ±0.30 
in distance m odulus. Th i s is a much larger value than 
that found by IKellv et all (|2010D . but the discrepancy is 
not significantly dramatic given our large uncertainties. 
Likewise, our v alue expressed m_dex/[0/H] is larger than 
that found bv iHowell etaLl (|2009l ) (-0.10 ± 0.07), but 
again the discrepancy is < 3er. For Cepheids, we find no 
significant correlation, but our measurement (0.23 ±0.22 
dex/mag or -0.92 ± 0.88 dex/[0/H] if the correlation 
arises from metallicity) can be compared to previous 
determinations. Given the large uncertainties, this re- 
sult agrees with previous values, whic h generally range 
betw een —0.1 and —0.4 dex/[0/H] (jRomaniello et all 
l2008h . 

3.5. The FM as Distance Estimator 

In Figure QT] we show the FM for galaxies that have 
their distances estimated using multiple techniques (we 
use the average for the adopted distance). The relation 
is satisfied by galaxies that range widely over luminos- 
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ity, morphology, and kinematics. It is tempting to invert 
the arguments above and use the FM not as a fiducial 
against which to test various distance estimators, but 
as a distance estimator itself. As stressed in the Intro- 
duction, these particular data are suboptimal for such a 
experiment, but we proceed nonetheless . 

Using the FM now for all of the galaxies in the sample, 
we estimate a distance to each, and show the relation- 
ship between recessional velocity and distance in Figure 
[T2l We define the FM zero point by setting the mean 
FM residual among estimators to zero, but define that 
mean in several different ways. First, we use the SBF 
calibration. Second, we calculate the weighted mean of 
the residuals in Table 2, excluding our calibration SBF 
sample. Third, we use the zero point from the SN la 
for pressure supported galaxies. The results differ by at 
most 4 km sec -1 Mpc -1 in the inferred Hq. We there- 
fore cite a systematic uncertainty of 4 km sec -1 Mpc -1 
in our Hq measurement below. We choose to plot the re- 
sults using the intermediate value of the FM zero point 
derived from the three methods (the mean of the resid- 
uals of the estimators). The inset shows the result of 
fitting the relationship for galaxies with v > 1500 km 
sec -1 , 3cr outliers excluded, binned by 10, and forcing 
the fit through the origin. We exclude the low v re- 
gion to minimize the effect of local flows. The resulting 
slope implies Hq — 78 ± 2 (random) ± 4 (systematic) 
km sec -1 Mpc -1 . The estimate of the systematic uncer- 
tainty does not include a variety of potential problems 
that are ignored here (modeling of bulk flows, internal 
extinction corrections, adjustment for potential biases in 
the galaxy sample, etc.). Therefore, while we treat this 
result cautiously, we present it to show 1) the level of 
precision possible with plausible sample sizes and 2) to 
illustrate how one could determine Hq in a way that uti- 
lizes as many (or as few) different distance estimators as 
desired. 

One of those sources of uncertainty, internal extinction, 
can be miti gated by going to the i nfrare d, as shown most 
recently bv lFreedman fc Madord (|2011[ ) for Spitzer wave- 
lengths where the scatter about the Tully-Fisher relation 
is reduced from 0.43 in the extinction-corrected B— band 
data to 0.31 in the non-extinction-corrected 3.6/zm data. 
We expect analogous gains in the FM for galaxies in 
which internal extinction is important. Depending on 
the fraction of galaxies in a sample for which extinction 
contributes significantly to the observed FM scatter, pro- 
portional gains in the determination of Hq will follow. 

4. CONCLUSIONS 

We demonstrate t hat the fundamental manifold (FM; 
iZaritsky et al.l 120081 ) provides a fiducial against which 
any set of distance estimators can be compared. Many, 
if not most, galaxies will only ever have distances mea- 
sured with one method, making it impossible to use 
them to test distance estimators directly. The FM fidu- 
cial, a benchmark for all distance estimators, constrains 
whether systematic errors exist among estimators and 
whether they are creeping into even a single estimator 
as a function of external parameters, such as metallic- 
ity Large sample size is critical for examining selection 
effects where, for example, one might want to split the 
sample in various ways (e.g., with environment, luminos- 
ity, morphological type, star formation rate). 




aooo 4000 

v [km sec -1 ] 



Fig. 12. — Hubble diagram using distances derived using the 
FM. Pressure supported systems are plotted as red circles and ro- 
tationally supported systems as blue squares. Inset shows data at 
v > 1500 km sec -1 binned so that each bin contains 10 galax- 
ies. The fitted line is constrained to go through the origin and 
corresponds to Hq = 78 ± 2 (random) ±4 (systematic) km sec 
Mpc -1 . 

Among the distance estimators we examine (SN la, 
Cepheids, surface brightness fluctuations, the luminosity 
of the tip of the red giant branch, circumnuclear masers, 
eclipsing binaries, RR Lyrae stars, and the planetary 
nebulae luminosity functions) , we find no statistically sig- 
nificant differences (> 2a), but caution that differences 
that are unacceptably large for state-of-the-art distance 
uncertainties (~ few %) are still allowed by our data. 
The lack of precision with which we could discriminate 
among methods is rooted in sample size, which in certain 
cases is still too small, and in the less-than-optimal data 
available (B band and heterogeneous). Both of these 
problems will be resolved with a feasible dedication of 
resources. 

We illustrate the use of the FM distances to 1) test the 
metallicity sensitivity of different estimators, confirming 
the SN la distance dependence on host galaxy metallic- 
ity and 2) provide an alternative calibration of Hq that 
replaces the classical ladder approach in the use of extra- 
galactic distance estimators with one that utilizes data 
over a wide range of distances simultaneously. 
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